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Abstract. The excitation of gravity waves by penetrative convective plumes is 
investigated using 2D direct simulations of compressible convection. The oscillation 
I field is measured by a new technique based on the projection of our simulation 

^■f^ . data onto the theoretical p-modes solutions of the associated linear eigenvalue prob- 

' lem. This allows us to determine both the excited modes and their corresponding 

, amplitudes accurately. 
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1. Introduction 



Two-dimensional simulations of compressible convection have shown 
that it is possible to excite internal gravity waves (IGW) in radiative 
zones of solar-type stars from the downward penetrating plumes (Hurl- 
burt et al- irnSFl hereafter HTM86; Hurlburt et al-. H^MI Kiraga and 
Jahn, I1995|) . However, detecting IGW with confidence is challenging 
given the stochastic nature of the excitation mechanism. 

We propose here a new detection method which allows us to measure 
rigorously both the spectrum and amplitude of excited 5-modes. This 
method is first applied to the 5- mode oscillations of an isothermal atmo- 
sphere and then, to IGW generated in a 2D-simulation of a convective 
zone embedded between two stable ones. 
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2. Detecting ^r-modes using the anelastic subspace 

In hydrosimulations, wave fields are commonly measured using two 
main methods: (i) the simplest one consists in recording the vertical 
velocity at a fixed point and then performing the Fourier transform 
of the sequence (see, e.g., HTM86); (ii) & more complicated method 
consists in taking two Fourier transforms, in space and time, of the 
vertical mass flux (Stein and Nordlund, I199U() . However, these two 
methods are not well adapted to detect IGW in our problem because 
the Fourier transforms are calculated over the whole simulation while 
IGW are stochastically excited by penetrating plumes. 

Our new detection method takes into account the random nature of 
this excitation. Indeed, it is based first, on projections of the simulated 
velocity field v{k,z,t) onto the anelastic eigenvectors ipkniz) as 

oo oo 
v{k, z,t) = ^ < V'fcn, V > tpkn{z) = ^ Akn{t)lpkn{z) , (1) 
n=0 n=0 

and second, on time-frequency diagrams of the complex coefficients 
Af:n{t)- As a consequence, the immediate spectrum (the set of fre- 
quencies uj) and amplitudes (defined as |Afcn|) of stochastically excited 
g-modes are reached and not only their "mean" values over the whole 
simulation. It is instructive to consider the simplest possible case, that 
is, the propagation of a single gravity mode with horizontal wavenum- 
ber k through the computation domain: applying eq. leads in this 
case to a projection coefficient Akn{t) oc eyip{iuJknt), where tOkn denotes 
the frequency of the anelastic eigenmode of degree k and order n. 



3. Results 

3.1. Oscillations of an isothermal atmosphere 

As a first test, we apply our method to detect IGW excited by an 
oscillating entropy bubble embedded in an isothermal atmosphere of 
depth d (see Fig. 1, left panel). In this case, the building of the anelastic 
subspace is simplified since we found analytic solutions for the eigen- 
frequencies and their associated eigenvectors (Dintrans et al. . 120021) . 

In Fig. 1 (right panel), we show the real part of the projection 
coefficient ^io(t), i.e. we projected the left panel velocity field onto 
the first anelastic eigenmode of the isothermal atmosphere at k = 1 
and 71 = 0. As expected, we found that Aio(t) behaves like exp(itJiot) 
(with LJio — 0.569cs/d, where Cg is the constant adiabatic sound speed) 
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Figure 1. Left: slide of the velocity field superimposed on the contours of entropy 
perturbations for a 2D-simulation of an oscillating entropy bubble embedded in 
an isothermal atmosphere. Right: time-evolution of the real part of its projection 
coefficient A\o (here t is in units of d/cs). 



whereas the mode amphtude |^io| decreases as exp(— with ty oc 
z/~^ {v being the constant kinematic viscosity of the simulation). 

3.2. 2D-SIMULATIONS OF PENETRATIVE CONVECTION 

Once our method vahdated, we study the excitation of IGW by over- 
shooting convection using high-resolution two-dimensional simulations 
of a three-layer polytropic model. That is, we solve the following equa- 
tions: 
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Dt 



div u, 



Du 

- = -(7- 
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where u denotes the velocity, e the internal energy, p the density, /C = 
K/cy where K is the radiative conductivity and the specific heat at 
constant volume, S the stress tensor and, finally, t(z) is a cooling time 
(see Brandenburg et al., ISM: for more details). 

Figure 2 (left panel) shows an example of such a simulation of a 
convective zone of depth d = 1 {0 < z < 1) embedded between two 
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Figure 2. Left: slide of the velocity field superimposed on the contours of entropy 
perturbations for a two-dimensional convection simulation. Two strong down- 
wards plumes penetrate into the bottom radiative zone and excite IGW. Right: 
time-evolution of the real part of the projection coefHcient A\_q (here t is in units of 



d/g). Thick lines emphasize IGW events detected using a time- frequency diagram. 



stable ones (—0.15 < z < and 1 < z < 3). IGW are excited in 
the bottom radiative zone by penetrating downward plumes and the 
evolution of the projection coefficient ^io(i) is now more chaotic (see 
right panel). However, by applying a time- frequency diagram on this 
sequence, we extracted three IGW events (with uiq ~ 0.251-^/ g/d), 
emphasized as thick lines in the figure. 
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